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ABSTRACT 



Aims. The aim of this study is to examine if the well-known chemical gradient in TMC-1 is reflected in the amount of rudimentary 
forms of carbon available in the gas-phase. As a tracer we use the CH radical which is supposed to be well correlated with carbon 
atoms and simple hydrocarbon ions. 

Methods. We observed the 9-cm A-doubling lines of CH along the dense filament of TMC- 1 . The CH column densities were compared 
with the total H2 column densities derived using the 2MASS NIR data and previously published SCUBA maps and with OH column 
densities derived using previous observations with Effelsberg. We also modelled the chemical evolution of TMC-1 adopting physical 
conditions typical of dark clouds using the UMIST Database for Astrochemistry gas-phase reaction network to aid the interpretation 
of the observed OH/CH abundance ratios. 

Results. The CH column density has a clear peak in the vicinity of the cyanopolyyne maximum of TMC-1. The fractional CH 
abundance relative to H2 increases steadily from the northwestern end of the filament where it lies around 1.0 10"**, to the southeast 
where it reaches a value of 2.0 10"^. The OH and CH column densities are well correlated, and we obtained OH/CH abundance ratios 
of ~ 16 - 20. These values are clearly larger than what has been measured recently in diffuse interstellar gas and is likely to be related 
to C to CO conversion at higher densities. The good correlation between CH and OH can be explained by similar production and 
destruction pathways. We suggest that the observed CH and OH abundance gradients are mainly due to enhanced abundances in a 
low-density envelope which becomes more prominent in the southeastern part and seems to continue beyond the dense filament. 
Conclusions. An extensive envelope probably signifies an early stage of dynamical evolution, and conforms with the detection of a 
large CH abundance in the southeastern part of the cloud. The implied presence of other simple forms of carbon in the gas phase 
provides a natural explanation for the observation of "early-type" molecules in this region. 

Key words. Astrochemistry - ISM: abundances - ISM: TMC-1 



1. Introduction 

The 9-cm A-doubling lines of the methylidyne radical CH 
have been commonly used to trace low density gas in diffuse 
clouds and at the edges of dense clouds. Observational surveys 
have demonstrated a linear correlation between the CH column 
density and visual extiction (and thereby also Af(H2)) in low- 
extinction molecular cloud s (e.g. iHjalmarson et a .ll977tlMattilal 
Il986t iMagnani & Onelld ll993l:TMagnani et al.l 120051 On the 
other hand, CH is likely to be critically dependent on the pres- 
ence of free carbon in the gas phase. The modelling results of 
iFlower et al. ] ( fTggI predict that CH (or CN) can trace atonuc 
carbon content even in the cloud interiors. While the fine struc- 
ture lines of CI become easily optically thick and dominated by 
the outer layers of the cloud, the ground state CH lines are opti- 
cally thin. Therefore one can hope to "see" through the cloud in 
these lines, even if the cloud edges may have a significant con- 



tribution, and the CH abundance is likely to vary as a function of 
density. 

We have made observations of the A-doubling transitions 
of CH in its ground state along the dense filament of Taurus 
Molecular Cloud-1 (TMC-1). With the aid of these observa- 
tions we examine the predicted correlation of atomic carbon and 
cyanopolyynes and other unsaturated carbon-chain molecules in 
TMC-1. 

According to most chemistry models adopted for TMC-1 the 
cyanopolyyne (CP) peak in its southeastern end should represent 
an early stage of chemical evolution and h ave a high concentra- 
tion of carbon in the form of C or C^Je.g. fHirahara et al.lll992t 
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* This publication is based on observations with the 100-m telescope 
of the Max-Planck-Institute fiir Radioastronomie (MPIfR) at Effelsberg. 



iHartquist et al.|[T99 ^: 'Pratap et al. |T997l). Alternative views have 
been presented by Hartquist et al.l (1 19961 ) who discuss the pos- 
sibility that the CP peak is highly depleted and represents a 
chemically more advanced stage than the NH3 peak, and by 
iMarkwick etaH (l2000[|200l who suggest that carbon chemistry 
at the CP peak has been revived by evaporation of carbon bearing 
mantle species (CH4, C2H2, C2H4) following grain-grain colli- 
sions induced by the ion-neutral slip produced by the passge of 
an Alfven wave. In this model, complex organic species are not 
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necessarily coincident with the most radimentary forms of car- 
bon. 

Previous CI 49 2 GH z fine structure line o bservations of 
ISchiUce et all (Il995h and iMaezawa etaP (1 19991) indicate high 
abundances neutral carbon in the direction of TMC-1. In their 
mapping of Heiles Cloud 2 (including TMC-1), Maezawa et al. 
estimate that the C I/CO abundance ratio towards the CP peak is 
about 0.1, and that the ratio increases to ~ 1 in the "CI peak" 
on the southern side of the cloud. These estimates are, however, 
uncertain because of the unknown optical thickness and excita- 
tion temperature of the 492 GHz CI line. The TMC-1 ridge is 
not visible on the C I map which is dominated by emission from 
more diffuse gas. 

In the present paper, we derive the fractional abundances of 
CH along an axis aligned with the dense filament of TMC-1 
and going through the CP peak. The molecular hydrogen col- 
umn densities, A^(H2), required for this purpose, are estimated 
by combining an extinction map derived from the 2MASS data 
with JCMT/ SCUBA dust contin uum maps at 850 and 450 //m 
presented in iNutter et"aD (l2008h . The CH observations are de- 
scribed in Sect. 2, and in Sect. 3 the CH and H2 column densi- 
ties are derived. In Sect. 4 we compare the new data with pre- 
vious observations, and discuss the chemical implications of the 
obtained CH abundances. Finally, in Sect. 5 we summarize our 
conclusions. 

2. Observations 

2.1. CH spectra 

The A-doubling line of CH in the '^H\i2,J =1/2 ground state 
{A - 9 cm) was observed in August 2003 using the 100-m ra- 
dio telescope, located in EfFelsberg, Germany, belonging to the 
Max-Planck-Institut fiir Radioastronomie. The acquired data in- 
cludes observations of the main spectral line 1-1 and the satellite 
lines 1-0 and 0-1. The frequency is 3335.481 MHz for the main 
fine, 3263.794 MHz and 3349.193 MHz for the lower and up- 
per satellites. The FWHM of the antenna at these frequencies is 
approximately 3'. 8. During the observations the system temper- 
ature varied approximately between 55 K and 60 K. 

The signal was fed into a 1024 channel autocorrelator which 
operated as a 3-band system, with the main line occupying 512 
channels with a bandwidth of 391 kHz. The satellite lines oc- 
cupied 256 channels each with a bandwidth of 195 kHz. With 
this setup the frequency resolution of each band is Av = 763 Hz, 
translating into a velocity resolution of Ai'n w 69 m/s for the 
main line. Cal ibration of the o bservations was done using flux 
densities from lOtt et al.l(ll994h of the well-known radio calibra- 
tors 3C84 and 3C123. The antenna temperature was observed to 
not be dependent on elevation as expected at these wavelengths. 

The observations were done using frequency switching. The 
origin of the observations (the 0,0 off'set) was set at q'j2ooo - 
04Hr"lZ7, (5j2ooo = 25°50'04" (o-Bigso = 04h38™08-^;7, (Jbi950 = 
2 5°44'20") which cor responds to the C'^O maximum in the map 
of lLanger et al.l(ll995h . The observations followed a line through 
this point at a 54° .5 tilt relative to the B1950 declinat ion axis. 
This li ne passed over the cyanopolyyne maximum C Olano et al.l 
(119881) . located near our position A.)c = +12') and the observa- 
tion poi nts partially overlapped with the previous OH observa- 
tions bv lHariu et al] (l2000t) . A total of 22 different off'sets were 
observed at the 3 wavelengths, bringing the total number of ac- 
quired spectra to 66. The offsets ranged from Aa - -8'. 13, 
A6 = +11'.39 to Aa = +16' .26, A5 = -22'.8 . See Figs.[Ta]and 
[lb] for a visual representation and Table [T]for a list of the off- 



sets. In addition to presenting the locations of the lines of sight 
as offsets along the right ascension and declination axes we will 
present them with Ax, which represents direct angular distance 
from the center of observations or (0,0). An example of reduced 
spectra is presented in Fig. |2] Due to human error, some of the 
initial observations were done using a slightly erroneous rest fre- 
quency. These observations had to be corrected manually back to 
the right rest frequency during the reduction phase. 

Furthermore, during our analysis we noticed that the CH 
emission peaks at an anomalous velocity of roughly 5.2 km/s. 
As no other observed molecules in TMC-1 are known to peak 
so prominently at this velocity and along these lines of sight, 
we initially considered the possiblity of CH tracing an unknown 
gas component in the line of sight of TMC-1. However, the fact 
that the obtained CH profiles are shifted with respect to OH 
aroused a suspicion that the CH velocities are incorrect as the 
two molecules are known to be well correlated in diffuse gas. 
We are grateful to the referee for pointing this out. 

After finding no explanation for the possible error in the 
spectral headers we looked for previous CH observations for ref- 
erence. The surroundings of the CP peak in TMC-1 have been 
previo usly observed in the CH main line by iMagnani & Onellol 
(Il995h using the NRAO 43 m telescope, and these spectra were 
kindly provided for us by Loris Magnani. gaussian fits to two 
of the Magnani & Onello positions closest to our observations 
provided velocity differences of ~0.5 km/s when compared to 
similar gaussian fits to our own observations. The peak veloci- 
ties of the Magnani & Onello positions were also similar to other 
molecules observed in TMC-1. The beam size of the 43 m tele- 
scope is 9' at this frequency and the velocity resolution of the 
spectra is 0.44 km/s. 

During the observing run we also did some test observations 
towards LI 642 (MBM20), wh ich has been previously observed 
in CH bv lSandell et al.l (Il98lh using the Onsala 25 m telescope. 
Our CH main line spectrum towards the (0,0) of Sandell et al. 
has a peak velocity of -0.22 + 0.05 km/s whereas the Onsala 
spectrum peaks at a slightly positive velocity: +0.2 ±0.1 km/s. 

After these comparisons we arrived to the conclusion that 
our CH lines have a spurious velocity shift which corresponds 
to about six autocorrelator channels. The effect of this shift is 
demonstrated in Fig.|3]which shows the CH and OH line profiles 
in various locations using both the original spectra as they come 
from the telescope, and spectra shifted by 6 channels (to an about 
0.4 km/s larger velocity). According to the EfFelsberg telescope 
team an error of this magnitude can have occurred in the con- 
version of the raw telescope data spectra to the Class reduction 
program format we use, in case an offset due to the quantization 
of the synthesizer frequencies in the IF section was not correctly 
registered. 

3. Results 

3.1. Calculation CH column densities 

Due to the very small Einstein A-coefficients of the transitions, 
the A-doubling lines of CH are optically thin. We can further- 
more assume that practically all CH molecules are in the ground 
rotational state and that the A-doubling energy levels are pop- 
ulated according to their statistical weights. The latter of these 
assumptions is valid because the separations between the levels 
are only ~ 0.16 K. 

Absorption measurements have shown that for CH 
the ground state A-doubling transitions are weak masers 
dRvdbeck et alJll975tlHialmarson et al.lll977h . In the case of the 




(a) (b) 

Fig. 1 : TMC-1 as seen by SCUBA at 850 (convolved to 30" in a and 228" in b). The map is reproduced from lNutter et al.l(l2008l ). 
who calls the filament "the Bull's Tail". The red diamonds represent the locations of our CH spectral observations (separated by 2'), 
and the black triangles m ark the locations of th e cyanopolyyne and ammonia maxima. The letters A-E and X mark the regions of 
the six dense cores as per lHirahara et al.l(ll992l ). with the asterisks marking the locations of the CCS maxima corresponding to the 
six cores. The two contour line levels correspond to intensities of 0. 19 and 0.32 Jy/beam on the 30" resolution map and 3.5 and 9.6 
Jy/beam on the 228" map. The filled red diamond on the image is located at ffjaooo = 04*'41™12';7, dmoo = 25°50'04", marking the 
origin of the CH observations performed by us. 



main line, the excitation temperature Tex 1 1 for CH is in the range 
~-60^-10K. 

Armed with this knowledge and the integrated intensity of 
the CH main line, J Tusdv [Kkms '], the total column density 
(in cm"^) of CH can be approximated with the formula 

jV(CH) ^ 2.82 10'^ \ r TMB(3335)dv (1) 

where rMB(3335) is the main-beam brightness temperature of 
the main A-doubling component (1 - 1), Tbg is the cosmic back- 
ground temperature and Texji is the excitation temperature of 
the 1-1 transition. 

In reality the assumption of statistical population distribution 
of the A-doubling energy levels is not precisely valid, which can 
be seen from the relative intensities of the three line components 
(Fig.|2j both satellites should be two times weaker than the main 
component). The total column density can be given as a function 
of the upper F=l column density Nu,i and the excitation temper- 
atures Tex,!!, Texjo and rex,oi: 

/ ''11 1 ''10 1 J±!__Jbl_\ 

A^(CH) = A^u.i 1 +■ e"«ii + -e''«.io + -e''cx,ii ''cx.oi 1 (2) 

We have used the observed relative integrated intensities of the 
line components, denoted by K^ y - J T^B.^dvl J TMB,ydv to 



derive the excitation temperatures Tex,io and Tgx.oi as functions 
T^ex.i I ■ The relative integrated intensities stay approximately con- 
stant throughout the observed filament. We obtain Kn^io ^ 2.30 
and /Til 01 ~ 1.44 by taking an arithmetic mean of their respec- 
tive area ratios from all the data. 

A value for Texji must be chosen for us to be able to get 
the CH column density. We c hoose Tex.ii = - 1 K w hich is 
similar to the value derived by iHialmarson et al.l (Il977b for the 
dark cloud LDN 1500 in front of 3C123. This choice implies 
that Texjo ~ -28.8 K and Tex.oi ~ ~3.5 K. Using the value of 
Tex,ii - -60 K from lGenzel et al.l (Il979h would result in an ap- 
proximately 23 % increase in the calculated CH column densi- 
ties. Using equation[T]for the column density calculations would 
have produced only a 0.3 % difference to these results. 

3.2. A closer look at the CH spectra 

The CH profiles are asymmetric in the sense that they peak 
close to the lower boundary of the radial velocity range of the 
detectable emission (see Fig. |2]i. It was noticed that at several 
offsets it is possible to fit two gaussian profiles to the observed 
spectra, especially for the main CH line: a narrow line compo- 
nent peaking between 5.4 and 5.7 km s"' , and a broader, weaker 
component roughly between 5.8 and 6.4 kms"'. In addition, a 
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Fig. 2: Some of the reduced observations in the vicinity of the 
cyanopolyyne maximum for each observed CH spectral line. 
The numbers on the left indicate the right ascension and dec- 
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Fig. 3: An example of overlapping CH (thick red) and OH (thin 



lination offsets in arcminutes from the center of the observations black) specti-a at several offsets before (left) and after (right) the 



at aBi950 = 04h38"'08^7, 631950 = 25°44'20" 



CH spectra have been shifted to a velocity 6 channels higher than 
in our reduced data. 



weak pair of wings can be discerned at some locations. The ve- 
locity distribution along the cloud axis is illustrated in Fig. |4] 
which shows the position-velocity diagram, and in Figs. ISaland 
|5b]presenting the results of the gaussian fits. 

We see from Fig.|5a]that the two gaussian components in the 
CH spectra can be resolved between offsets -4' . . . + 18' along 
our axis of observation. Looking at the points in the SCUBA 
maps we see that the offset Ax = H-18' lies near the southeastern 
tip of the filament, whereas Ax = -4' is already somewhat off 
from the bright part on its northwestern side. 

Of the aforementioned two gaussian components the more 
narrow component starts off at about 5.4 km/s and then suddenly 
jumps to Vlsr ~ 5.7 km/s at A.y = +8' and again to about 5.8 
km/s at Ax = +20'. The peak of the broader component persists 
at LSR velocities roughly between 6.0 and 6.1 km/s. As can be 
seem from Fig. |5b] the low-velocity components also seems to 
become significantly narrower at Ax = +10', halving its FWHM. 
This halved line width is equivalent to thermal broadening (as 
calculated from dust temperatures in Table [1} at these locations, 
implying that the CH emission of that line is emitted from an 
area of the cloud where no other line broadening effects, e.g. 
turbulence, are present in any significant part. We note that the 
narrowing of masing lines is significant at large negative optical 
thicknesses, -t » 1 (e.g. Spitzer 1978, Ch. 3.3.), but for CH 
lines with |t| < 0. 1 the effect is negligible. 



At other offsets the width of the sharper line can not be en- 
tirely accounted for by thermal broadening and the wider line 
seems to have a significant component of non-thermal broaden- 
ing. Overall the non-thermal width of the sharp line varies be- 
tween about m/s and 400 m/s while the non-thermal width of 
the wide line persists between 1 km/s and 1.4 km/s. 

Due to the peculiar narrow blueshifted component, we de- 
cided to examine the calculated CH column densities not only in 
the full velocity range of detectable emission (4.5-6.5 km/s) but 
also with velocity ranges roughly corresponding to the two sep- 
arate line components. We calculated A^(CH) for the Vlsr range 
4.5 - 5.4 km/s for the narrow blueshifted component and for 
the range 5.4-6.5 km/s for the wider component. All three be- 
haviours of A^(CH) along the filament are presented in Fig. |6]and 
the total column densities of CH are presented in number form 
in the fourth column of Table [1] 

In the figure we see that the narrow component seems to be 
the major contributing factor to the slow increase of the total 
A^(CH) starting at around Ax = -6', and its likewise gradual 
decrease at around Ax = +20'. The A^(CH) contribution of the 
broad component on the other hand seems to stay almost con- 
stant throughout the entire filament, increasing only very slightly 
as we go from the northwestern end of the filament to the south- 
eastern end. It does, however, have a distinctive "bump" in the 
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Fig. 4: Position-velocity diagram of the main CH spectral line. 
The lowest contour line corresponds to 0.27 K, increasing by 
about 0.05K per contour 
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Fig. 5: (a) Centroid positions and (b) widths of the two gaussian 
fits made to the main CH line at different offsets. The triangles 
represent the blueshifted component of the two fits. 



Fig. 6: Column densities for CH along the observed lines of sight 
in TMC-1. The black diamonds represent column densities cal- 
culated from the full velocity range and the blue triangles and 
red squares represent column densities calculated from velocity 
ranges corresponding to the narrow and wide gaussian compo- 
nents respectively. 



vicinity of the cyanopolyyne maximum and it seems to be the 
major contribution to the similar bump seen in the total A^(CH) 
at that region. 

This bump is obviously also the maximum of the total col- 
umn density in the entire region we measured and it is located 
approximately at Ax = 14', roughly between Hirahara's cores D 
and E. 



3.3. H2 column densities 



The NICER method dLombardi & AlvesI I200TI) was appUed to 
2MASS data of the filament region of TMC-1 in order to acquire 
visual extinction values for the points overlapping with our CH 
observations, with the angular resolution element corresponding 
to the FWHM available in the radio observations. Then, by ap- 
plying the E(J - K)/N(U2) - 5. 4 10^' cm'^mag ' rela tion from 
iHarjunpaa & Mattilal (Il996h and lCardelli et al.l (Il989h . it would 
normally be possible to get column densities for molecular hy- 
drogen from the extinction data alone. However, as the dense 
cores along the filament have high extinction values, acquiring 
reliable values for Ay around them by using the NICER-method 
alone is unreliable due to very low star densities close to and 
inside the cores. 

As an alternative means of obtaining H2 column densities 
we used submillimetre dust continuum map s. The 850 and 450 
IJ.m SCUBA maps of the TMC-1 region from lNutter et al.l(l2008h 
were kindly provided for us by David Nutter After convolving 
them to the FWHM of our 9 cm observations (228") these were 
used to estimate A^(H2) in the positions observed in CH. The 
scan-mapping mode used in the SCUBA observations involves 
chopping by a maximum throw of 68". The continuum maps are 
not sensitive to spatial scales larger than that. In order to com- 
pare the dust maps with the extinction and CH line data, we need 
to convolve the 850 and 450 mum SCUBA maps to the FWHM 
of 228" of our 9 cm observations. We believe that in the vicinity 
of the TMC-1 ridge where the dust emission is likely to be dom- 
inated by compact cores this convolution results in reasonable 
"smoothed" flux densities. Furthermore, we attempt to correct 
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Table 1: Column densities of different molecules, visual extinctions and dust temperatures along the TMC-1 filament. The velocity 
interval of the integrated main beam temperature used for calculating A^(CH) and A^(OH) is 4.5 km/s . . .7.5 km/s. Ao- and A<5 
represent right ascension and declination offsets and Ax represents direct angular distance from the center of the observations. 



Ax ['] 


Aa ['] 


A5['] 


A'(CH) [10" cm-'] 


W(OH) cm-^] 


Av [mag] 




[K] 


A'CHz) [10- 'cm--] 


+28.00 


+ 16.26 


-22.80 


1.11 ±0.08 


2.27 ±0.12 


4.82 ±0.18 




- 


4.53 ±0.17 


+26.00 


+ 15.10 


-21.17 


1.12 ±0.08 


- 


4.90 ±0.19 




- 


4.61 ±0.18 


+24.00 


+ 13.94 


-19.54 


1.26 ±0.07 


2.12 ±0.07 


5.12 ±0.18 




- 


4.81 ±0.17 


+22.00 


+ 12.78 


-17.91 


1.28 ±0.08 


- 


5.89 ± 0.24 


12.61 


±0.37 


5.54 ±0.23 


+20.00 


+ 11.61 


-16.28 


1.38 ±0.08 


2.50 ±0.11 


7.55 ± 0.37 


14.69 


±0.38 


7.10 ±0.35 


+ 18.00 


+ 10.45 


-14.65 


1.34 ±0.07 


- 


8.46 ± 0.52 


15.49 


±0.37 


7.95 ± 0.49 


+ 16.00 


+9.29 


-13.03 


1.61 ±0.07 


2.74 ± 0.08 


9.74 ± 0.64 


14.12 


±0.30 


9.15 ±0.60 


+ 14.00 


+8.13 


-11.40 


1.66 ±0.08 


- 


12.20 ±0.89 


13.04 


±0.26 


11.47 ±0.84 


+ 12.00 


+6.97 


-9.77 


1.61 ±0.07 


2.69 ± 0.05 


11.00 ±0.82 


12.87 


±0.26 


10.34 ± 0.77 


+ 10.00 


+5.81 


-8.14 


1.48 ±0.07 


- 


12.51 ±0.92 


12.77 


±0.26 


11.75 ±0.87 


+8.00 


+4.65 


-6.51 


1.37 ±0.08 


2.68 ± 0.07 


14.04 ± 1.05 


12.14 


±0.25 


13.20 ± 0.98 


+6.00 


+3.48 


-4.88 


1.37 ±0.08 


- 


27.01 ± 7.39 


11.52 


±0.25 


13.11 ±0.55 


+4.00 


+2.32 


-3.26 


1.37 ±0.07 


2.54 ± 0.08 


120.92 ±73.83 


11.24 


±0.24 


13.42 ± 0.58 


+2.00 


+ 1.16 


-1.63 


1.40 ±0.07 




10.48 ±53.10 


11.66 


±0.26 


12.05 ± 0.53 


0.00 


0.00 


0.00 


1.24 ±0.08 


2.54 ± 0.06 


11.58 ±0.82 


10.95 


±0.28 


10.89 ± 0.78 


-2.00 


-1.16 


+ 1.63 


1.18 ±0.08 




9.72 ± 0.62 


8.73: 


+ 0.28 


9.13 ±0.58 


-4.00 


-2.32 


+3.26 
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for the arbitrary background level of the SCUBA maps by cor- 
relating the intensities with the extinction values as desribed be- 
low. The 850 jum SCUBA maps around the Bull's Tail smoothed 
to 30" and 228" resolutions are shown in Figs. [Ta] and [lb] re- 
spectively. 

The H2 column density is related to the intensity of the dust 
emission by 



By(Td) ■ 2.^mH ■ Ky ■ R/ 

where By{TS) is the Planck function, is the temperature of the 
dust, mil is the mass of a hydrogen atom, Ky is the dust cross- 
section per unit mass of dust and is the dust-to-gas ratio. For 
Ky wdRd we used the values 0.13 cm^g-' (at 850 //m) and 1/100, 
respectively. 

As bolometer maps usually do, the SCUBA maps con- 
tain negative artefacts, i.e. hollows, around bright emission re- 
gions. This causes an uncertainty to the absolute level of the 
intensity /,,, and suggests a need for performing a bias correc- 
tion. We combined the NICER-generated extinction map with 
the SCUBA maps to examine the relation between extinction 
and submm intensit y much in the same way as was done by 
iBianchi et al.l (l2003h . This examination can provide us with two 
things: A bias correction for the SCUBA maps, and the average 
temperature of the dust along the filament. 

A linear fit of the form 

if'^Sy+JyAy. (4) 

was performed to the pixel by pixel correlations between the ob- 
served intensities at 850 or 450 and the visual extinction. 
The fit was done using pixels with Ay < 12, because above this 
level the extinctions are highly uncertain (see above), and be- 
cause the dust temperature, Td, is expected to decrease towards 
the most obscure region s, causing a curvature in the correlation 
(see lBianchi et al.ll2003h . 



The coefficient 6y should be zero in a properly biased map. 
The intensity of the dust emission is given by 

/f^' = By(Ti){\ - e-^'O * By{Ti)Ty , (5) 

where the optical thickness, Ty, of the source of emission can 
be written as Ty - l,KyRd, where Z is the surface density of the 
cloud. It can be assumed that subtracting the fitting parameter 
6y from the 850 and 450 /im intensity maps should result in a 
better- than-no thing correction to the bias. 

The dust temperature was obtained from the ratio of the y 
coefficients at the two frequencies: 

Jyl _ Kyi Byl(Td) _ B yl(T d) 

JVI Ky2By2(Td) \V2/ By2(Td) ' 

In the latter form, we have made use of the assumption that at 
(sub)millimetre wavelengths Ky cc w here the emissivity in- 
dex y6 - 1 .95 between 850 and 450 //m (lOssenkopf & Hennind 
IT994I) . 

As the dust temperature, and possibly also the ratio dust 
opacities in the submm and in the visual, Ky/Ky, can be expected 
to be different in the filament when compared to several other 
parts of the cloud, the examination of the linear relation between 
observed intensities and extinctions were done only on pixels 
in a rectangular area passing along the filament with a width of 
about 400". Linear fits were performed in these areas for com- 
paring both the 450 /^m and 850 //m maps to the visual extinc- 
tions as calculated with the NICER method. These comparisons 
and the linear fits performed on them are presented in Figs. |7a] 
and|7bl With the help of these fits we corrected the bias on the in- 
tensity maps, and calculated the average dust temperature of the 
examined area to be = (1 1.0 + 2.7) K. After this, with the help 
of Eq. [3] it was possible to use the 850 /im (because it seemed 
more accurate) SCUBA map to calculate the H2 column densi- 
ties on the points overlapping with our CH observations (three 
points in the southeast were not covered by the SCUBA maps). 
The resulting A^(H2) estimates are plotted in Fig. |8]as a dashed 
Une. 
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The column densities calculated in this fashion agree rather 
well with those estimated from the Ay map on the northwestern 
side of the map {Ax < 0', see Fig. [8]), whereas the agreement is 
not good in the southeast. However, it is conceivable that the Ay 
map reflects the true H2 column density distribution except for 
the regions of the highest obscuration, in particular it should do 
so at both ends of the strip observed in CH. 

The discrepancy between the N(H2) estimates from the 850 
fim emission presented above and those from Ay can be caused 
by the fact that assumption of a constant dust temperature is not 
valid. We examined therefore the possible variations of the dust 
temperature along the filament by using the 'bias-corrected' 450 
jum and 850 /im intensities in the formula 



850/xm, FWHIv1=228" (masked) 



\V2 



3.0+/3 hvJkT, 



1 



(7) 



This method can be applied with a reasonable accuracy to the 
positions within the dense filament where both 450 and 850 yum 
signals are strong. The resulting dust temperatures and H2 col- 
umn densities are presented in Table [1] and as a dotted line in 
Fig. [8] It is worth noting that the that, in the envelope, the dust 
temperature calculated in this fashion is similar to the = 12 K 
derived by Nutter et al. for the shoulder component. This is also 
quit e close to kinetic g as temperatures observed in the cloud by 
e.g- lPratapetaH (119971) . 

Since the extinction method seems reliable everywhere ex- 
cept at three points (Ax -2' . . .6') around the calculated extinc- 
tion value maximum and because the dust temperature seems 
to be variable along the filament, we have opted to use the 
extinction-derived A^(H2) values for all except the three afore- 
mentioned points. For these three points we use instead the 
A^(H2) values calculated by using dust temperatures acquired 
from Eq. |7] These column densities are listed in the sixth col- 
umn of Table [U 
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3.4. CH abundance 

When comparing the column densities of H2 and CH in an x,y- 
plot, the slope of a linear fit performed on the data points should 
tell the average abundance of CH in the observed area. This ap- 
proach is however not very interesting if the two column densi- 
ties do not correlate very well. No strong correlation (p w 0.53) 
between the two column densities was noticed, suggesting that 
the fractional abundance of CH is not very constant throughout 
the region. The abundance of CH, or X(CH) = A^(CH)/A?(H2), 
along the observed axis is presented in Fig.|9] This figure shows, 
beginning from Ax = +8', a steady rise in X(CH) from about 
1.0 10 up to ~ 2.2 10"^ towards the southeastern end of the 
filament. Our result at the cyanopolyyne peak matches the value 
quoted in Ohishi et al. ( 1992). Examining the abundance of CH 
calculated by using A^(CH) from the narrow and wide line com- 
ponents is not very useful because we can not evaluate A^(H2) 
from the two different velocity intervals separately and thus we 
refrained from doing so. 

3.5. Correlation with OH 

We had access to previous OH spectral line data acquired by 
iHariu et aH (I2OOOI) with partially overlapping points to our CH 
observations and used the formula from the same paper to ac- 
quire the column densities for OH: 

iV(OH) = 4.04 10'4— — 1— . — r rMB(1665)flfv (8) 
1 - Tbg/rex rjc J 



(b) 

Fig. 7: A linear fit to the visual extinctions and (a) 850 yum or (b) 
450 yum intensities around the filament region of TMC-1. The 
resolution of the intensities and the extinctions both correspond 
toabeamFWHM of 228". 

with [Kkms"'] as the unit for the line area and [cm"-] as the 
unit for the resulting column density. The excitation temperature 
and main beam efficiency used for this was Tex = 10 K and 
rjc - 30 %. The velocity interval (4.5 km/s . . . 6.5 km/s) of the 
integrated line area was chosen to be identical to the one used 
for the CH column density calculations. The resulting column 
densities are listed in the fifth column of Table [T] and they are 
plotted in Fig.fTOl 

As can be seen from Table [1] there are twice as many CH 
observations as there are OH observations along the strip we are 
investigating. To make use of all the CH observations when com- 
paring them to OH, we have used a simple but sufficient averag- 
ing scheme where, on each overlapping CH and OH observation, 
we take a weighted average of CH and its adjacent point on each 
side and give the adjacent points a weight of one half. The over- 
lapping observations are plotted to Fig.[TT]along with a linear fit 
made between the points. A strong correlation (/?ch,oh ~ 0.90) 
is noticed between the column densities and the column densi- 
ties follow the relation 

A^(OH)[cm-'] = (11.2 + 2.1) 10'^ + (10.1 ± 1.6)-A^(CH) . (9) 
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Fig, 8: Column densities for H2 along the TMC-1 filament. The 
dashed line represents molecular hydrogen column densities cal- 
culated with Eq,|3] while assuming a constant dust temperature 
of r^; » 11.0 K, The diamonds represent A^(H2) values derived 
from a NICER extinction map as described in the beginning of 
section 13.31 The triangles and the dotted line represent A^(H2) 
values calculated by using variable dust temperatures obtained 
from Eq. [T] 



X(CH) 
X(OH) - 
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Fig. 9: The abundances of CH and OH along the observed lines 
of sight in TMC-1. 



The fractional abundance of OH along the TMC-1 ridge 
is plotted together with that of CH in Fig. |9] As the CH col- 
umn density lies in the range ~ 0.9 - 1.6 10'"^ cm"^ the de- 
rived OH/CH abundance ratio along the ridge, A^(OH)/A?(CH)~ 
16 - 20, is clearly larger than tha t meas ured recently in diffuse 
interstellar gas bv IWeselak etaP (l2009l) : N{OH)/N{CH)~ 2.6, 
which is close to the atomic O/C ratio in the local ISM (e.g. 
iDulev & Wiliiarnsl[T98l I Wilson & Roodl[T99l . 



4. Discussion 

4.1. The chemistry of CH and its relation to OH 



Pure gas-phase chemistry models (e.g. IVialal[T986l) predict that 
the CH abundance, X{CH), decreases from ~ 10"^ at the outer 
boundaries of molecular clouds (Ay < 1) to < 10"^ towards the 



Fig. 10: Column densities for OH and CH along the TMC-1 fil- 
ament. The arrows point out the approximate lo cation of the 
cyanopolyyne maximum as per lOlano et al l (ll988h . 




Fig. 11: Correlation plot of the column densities of CH versus 
OH in spots where observations of them overlap. 

cloud interiors. This change should imply saturation in the cor- 
relation of A^(CH) with A^(H2) towards dense molecular clouds, 
and ind eed this kind of tendency can be traced in observation al 
results (lHialmarso netaLlll977t lMattilalll986l: lOin et al.ll2010l) . 

The production of CH depends on the availability of C^ or C. 
C^ is the dominant form of carbon in the outer parts of a cloud 
exposed to UV radiation. In the inner parts, C^ formation is hin- 
dered by the lack of UV photons and if formed is rapidly neu- 
tralized or incorporated into neutral molecular species via ion- 
molecule and dissociative recombination reactions. When the el- 
emental O/C ratio is larger than unity, most carbon is consumed 
by the formation of gas-phase CO, which eventually freezes onto 
grain surfaces at sufficiently low temperatures, < 20 K. 

The principal source of CH in dark clouds is likely to be 
the dissociative recombination of hydrocarbon ions, in particular 
CH+ and CH+, e.g.: 



CHt +, 



CH + Hz 



(10) 



The sequence leading to these ions is initiated either by the ra- 
diative association of C^ and H2: 



C^+Hz 



CH^ 



• hv 



(11) 
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or by the proton exchange reaction between neutral C and : 



C + CH+ + H2 



(12) 



(e.g. iBlack et alj [19751 and references therein). The sequence 
continues to larger hydrocarbon ions through reactions with H2. 
CH is thus expected to correlate with and C I, and because 
these species are necessary to the production of complex carbon 
compounds, it should correlate also with heavier hydrocarbons 
and cyanopolyynes. Furthermore, a close relationship between 
CH and long carbon chains or cyclic molecules may originate in 
direct CH insertion into an unsaturated hydrocarbon (followed 
by H elimi nation) which is like ly to proceed with a small or no 
barrier (e.g. lSoorkia et al.ll2010[ and references therein). 

The destruction of CH is thought to occur primarily via the 
following neutral exchange reactions: 



CH + 
CH + N 



CO + H 
CN + H . 



(13) 



In contrast with pure gas-phase chemistry, models in- 
cluding accretion of molecules onto grain surfaces predict a 
smooth distribution for the CH abundance in dense c louds. 
iHartquist & WiUiamsl (1 19891) and iRawlings et al] ( Il992h have 
pointed out that the CH abundance may remain almost un- 
changed in the dense interiors of dark clouds where heavy 
molecules are likely to be depleted. As long as some gas-phase 
CO is available, carbon ions are supplied by reaction with He^ 
(which is continually produced by cosmic rays): 



He+ + CO ^ C+ + O + He 



(14) 



The reduction of CO is partly compensated by the diminishing 
destruction of CH by O or N. Also an increase in the HJ 
abundance associated with the depletion can have a balancing 
effect on the CH abundance. Furthermore, the adsorption energy 
of CH is relatively low and it is efficiently return ed to the 
gas-phase via thermal desorption (lAikawa et al.lll997l) . 

In dark clouds, the reactions leading to the formation and de- 
struction of OH are probably analogous to those forCH. It is thus 
expected that OH is mainly produced by dissociative recombina- 
tion of H30^ and destroyed in neutral exchange reactions with 
C, N, and O atoms. In order to identify the most important re- 
actions and to study temporal changes in the OH/CH abundance 
ratio, we ran some chemical models which, although relatively 
simple, are appropriate for dark clouds such as TMC-1. 

We modeled the chemistry of TMC-I by treating the source 
as a homogeneous, isotropic cloud with the constant physical 
parameters, «(H2) - 10"* cm"-', T = 10 K (for both gas and 
dust) and Ay = 10 mag, o-g «; 0.03 jum, using a cosmic-ray 
ionisation rate of 1.3 10 '^s ' and performing a pseudo-time- 
dependent calculation in which the chemistry evolves from ini- 
tial abundances that are atomic, excepting molecular hydrogen. 
We a dopt the low-metallic ity oxygen-rich elem ental abundances 
from lGraedel et al.l(ll982l) a nd lLee et al.l(ll998l) as hsted in Table 
8 of Woodall et al ] (I2007h . Thus, our initial elemental abun- 
dances with respect to H nuclei number density for C and O are 
7.3 10"^ and 1.76 10""*, respectively, giving a C/O ratio of 0.4. 
Note, under dark cloud conditions, hydrogen exists primarily in 
molecular form so that the H nuclei number density is twice 
that of molecular hydrogen. Our gas-phase chemical network 
is the latest release of the UMIST Database for Astrochemistry 
or "Ra te06" ( http://www.udfa.net| l as described in I Woodall et al.l 
(I2007h . We ran one model using gas-phase chemistry only and 



< 

ra 
c 
o 



10"- 
10-' 
lO"- 
10-' 
10"' 
10-' 
10-^ 
10-" 
IQ-" 



Gas-phase Model 








1 








CH 

OH 

C 

CO 









10'' 



10"* 



10"" 
Time (yr) 

(a) 



10" 



10' 



10-- 



10-' 



s 10-- 



< 

c 
o 



Gas-grain Model 










CH 1 

OH 1 

c 1 

CO 11 



10-' 
10-' 
10-' 

u 

to 

10-- 

10-" 

10-^- 

10'' 10"* 10' 10° 10' 

Time (yr) 
(b) 

Fig. 12: Plots showing the fractional abundances (with respect to 
H2 number density) of CH and OH as functions of time, calcu- 
lated using a gas-phase chemistry only model (a) and a gas-grain 
model (b) for a TMC-1 -like cloud (see Sect. 14. lb . Also shown in 
this Figure are the fractional C and CO abundances. The shaded 
regions represent the observed range of CH and OH fractional 
abundances derived by us along the TMC-1 ridge. 



another allowing for gas-grain interactions, namely, the accre- 
tion of gas-phase species onto dust grains and the thermal des- 
orption of species from the grain surfaces. For our second model, 
we assume the dust grains are spherical particles with a radius of 
lO"*" cm and that the dust is well mixed with the gas possessing 
a constant fractional abundance of 2.2 10"'^ with respect to H 
nuclei density. Our accretion and thermal desorption rate s were 
calcul ated according to the theory outlined in Hasega wa et al.l 
assuming a sticking coefficient of unity. 
The fractional abundances (with respect to H2 number den- 
sity) of CH and OH are plotted in Figs. ll2a land ll2bl as functions 
of time as predicted by the pure gas-phase (solid lines) and gas- 
grain (dot-dashed line) models. We see that the two models agree 
rather well until an age of approximately 210^ years when the 
freeze out of molecules starts to take effect. In both models, the 
CH and OH fractional abundances show an anticorrelation un- 
til a time of about 5 10^ years, with CH continually declining 
and OH continually increasing. The OH abundance exceeds that 
of CH at an early time before 10"* years and remains thus for 
the duration of the simulation, except until very late stages in 
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the gas-grain model. CO passes C in abundance in these models 
around 6-7 10"* years. 

Soon after 10"^ years CH and OH plateau temporarily (until 
~ 2 10'' year). The relative constancy of the OH/CH abundance 
in this time interval can be explained by very similar gas-phase 
production and destruction mechanisms. 

According to the present model, the dissociative recombina- 
tion of HiO^ is responsible for 81 % of the OH production at a 
cloud age of lO"* years and for 89 % at 10^ years. It is only at 
the relatively young cloud age of 10^ years that dissociative re- 
combination contributes less than 20 %. At this stage, the model 
predicts the reaction 

O + HCO ^ CO + OH (15) 

to be a major process, making up 76 % of OH formation. 
Likewise, the dissociative recombination of CHj is accountable 
for 72 % of CH production at 10^ years and for 69 % at 10^^ 
years, yet contributes only 23 % at 10^ years. The main produc- 
tion pathway at this age is 

CH2+H^CH + H2. (16) 

Since CH2 is also likely produced through dissociative recom- 
bination of CH3 , this "detour" does not significantly change the 
overall picture. All these precursors are formed by successive 
hydrogenation of CH+ and CHj ions that stem from reactions 
of C^ with H2 or C with H^. Also, HjO^, the main precursor 
of OH, has a similar genesis through iterative hydrogenation of 
OH^ formed by protonation of O atoms by H^ . 

Furthermore, the destruction pathways of CH and OH are 
very similar. At a cloud age of 10^ years, reactions with C, N 
and O, 

OH-hC ^ CO + H(31%) (17) 
OH + N ^ NO-hH(19%) (18) 
OH + O ^ 02H-H(49%) (19) 

are exclusively destroying OH (the figures in brackets denote 
the relative importance of each reaction). CH is consumed by 
the following analogous processes: 

CH + C ^ C2+H(14%) (20) 
CHh-N ^ CNh-H(21%) (21) 
CHh-O ^ COh-H(60%) (22) 

A small proportion is, interestingly, due to the associative ion- 
ization (reverse dissociative recombination) of CH with O 

CH H- O ^ HCO+ + e" (23) 

which is almos t thermoneutral (alth ough if one calculates the 
Gibbs energies dHamberg et al.ll20Tol) . it becomes endoergic and 
probably should be removed from the model). At earlier stages 
(approximately 10^ years), reactions with C^, 

C+ -H CH ^ C^ -H H (42%) (24) 
C+ -H CH ^ CH+ + C (42%) (25) 

are more important, because of the higher abundances of C^ at 
this time. 

The pure gas-phase and gas-grain models start showing large 
deviations from each other at times later than lO*" years. The 
fractional abundance of CH declines quickly as C is converted 
into CO, and in the pure gas-phase model it diminishes contin- 
ually until a steady-state fractional abundance of » 3 10"" is 



reached at a time of about 2 10^ years. In the gas-grain model 
around a time of 10^ years, the revival of OH and CH is due to 
the depletion of atomic species (C, N and O) as their abundance 
drops rapidly after a time of about 5 10^ years thus limiting their 
effectiveness at destroying the radicals. Beyond 10^ years, stable 
species (e.g. CO, N2, O2 and H2O) become locked onto the grain 
surfaces thus depleting the gas-phase of all atoms and molecules. 
In the present model, which includes thermal desorption from 
grains only, both CH and OH abundances decrease rapidly af- 
ter about 2 10^ years. This behaviour is different from models 
which include non-thermal desorption mechanisms (cosmic ray 
heat ing, surface p hotodissociation, exothermic surface reactions; 
e.g. iGarrod et al.l (120Q7)) which predict detectable amounts of 
CH and OH at later times. 

4.2. Distributions of CH and OH in TMC-1 

When looking at our data we see in Fig.|6]a clear peak in the CH 
column density in the vicinity of the cyanopolyyne maximum of 
TMC-1, about 3' on its southeastern side (offset Ax = +14'). The 
bump in A^(CH) begins at the offset Ax = +8' which corresponds 
to the point where our strip enters the region of greatest intensity 
on the smoothed SCUBA map (Fig.[Tbl), and ends at Ax = +18' 
where we exit the brightest parts. 

The diagram in Fig.|6]shows also CH column densities in two 
velocity windows, and it is evident that the gas in the velocity 
range 5.8 - 7.5 kms"' is mainly responsible for the peak. Most 
molecular lines observed towards this source, including those 
of complex carbon compounds, originate in this velocity range 
(e.g. Pratap et al. 1997), and it seems likely that the local CH and 
cyanopolyyne peaks are also nearly spatially coincident. Besides 
the peak at Ax - +14' the A^(CH) in the Vlsr range 5.8 - 7.5 
kms"' shows a tendency to increase slightly from northwest to 
southeast, all the way until the end of the strip where the dust 
column density and A^(CH) in the low-velocity range, 4.5 - 5.8 
kms"' have decreased. 

With respect to the fractional CH abundance, X(CU), the CP 
peak is not special. From Fig.|9]we see that X(Cli) lies some- 
where around 1.0 10 in the northwestern part of the filament 
until the offset Ax = +8', and reaches values of ~ 2.0 10"** to- 
wards the southeastern end. This is the only major trend in the 
variation of X{CH) and it does not seem to vary in e.g. the vicin- 
ity of the high-extinction cores. 

It is evident from Figures l9lfTn that the CH abundance is well 
correlated with that of OH. The fractional OH abundance in- 
creases from about 2.5 10"^ in the northwest to about 5.0 10"^ in 
the southeast. Acccording to the chemical models discussed in 
the previous Section, a good correlation between OH and CH is 
expected at youthful cloud ages around 10"* - 10^ years because 
of the very similar gas-phase production and destruction mech- 
anisms. During this period, the OH/CH abundance ratio in the 
models increases from about 3 to about 10. It may be possible to 
reproduce the observed OH/CH ratio by adjusting the elemen- 
tal O/C ratio and/or the physical parameters of the model. We 
note that during the model ages between 10"* and 10^ years the 
fractional C and CO abundances are within one order of magni- 
tu de from each other, whic h roughly agrees with the estimates 
of lSchilkeetaP (Il995h and lMaezawa et all (Il999l) . 

On the other hand, the gas-grain model predictions repro- 
duce the observed large abundances of CH and OH at later times, 
around 10* years. Furthermore, the gas-grain model can qualita- 
tively explain the observed spatial variations of the CH and OH 
abundances along the TMC-1 filament if they are interpreted in 
terms of evolution with the chemical age advancing from the 
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southeast to the northwest. After the peak in OH in the gas-grain 
model (~ 1.5-2 10^ years), both CH and OH fractional abun- 
dances decrease with increasing age, and it is only at this point in 
the present calculations that this particular trend is seen. Running 
models including a more comprehensive treatment of desorption 
may shift this period but most likely it would still be associated 
with an advanced stage where accretion onto grain surfaces plays 
a dominant role. 

In view of the rather coarse angular resolution available in 
the present observations, and the fact that both CH and OH are 
known to have large abundances in diffuse gas, it is, however, un- 
likely that observed abundance gradients predominantly reflect 
conditions in depleted cores, which occupy a relatively small 
fraction of the gas encompassed by the telescope beam. A more 
probable explanation is that when moving from the northwestern 
end of the filament to the southeast we see an increasingly thick 
layer of tenuous gas, which would also agree with the increasing 
abundance of neutral carbo n towar ds southeast suggested by the 
CI map of iMaezawa et al.l (Il999l) . This explanation is further- 
more supported by our chemical models where, if we decrease 
the gas density, we can see the abundances of both CH and OH 
increase. In this situation the density dependence overcomes the 
temporal anticorrelation. This tendency is probably due to an en- 
hanced production by dissociative recombination resulting from 
the higher fractional abundance of electrons in low-density gas. 

[Hirahara et al. ( 1992) and Pratap et al. ( 1997) estimated by 
modelling the excitation of molecular lines that the gas volume 
density is somewhat higher near the ammonia peak in the north- 
western part of the TMC-1 filament as compa red with the CP 
peak. The SCUBA maps of iNutter et"aD (l2008h , reproduced in 
the present work, show that the CP peak lies in the area where the 
dust and gas column density reaches its maximum. This separa- 
tion of volume and column density peaks implies that the south- 
eastern part of the filament possesses a more extensive envelope 
of low-density gas than the northwestern part, and is readily an 
indication that the stages of dynamical evolution are different at 
the two ends of the filament. 

The large abundances of CH combined with the presence of 
the SCUBA dust intensity maximum in the southeastern parts 
is consistent with the presence of the low-density envolope. The 
fractional abundance reached at the southeastern end of the strip, 
X(CH) ~ 2.0 10 is slightly below half when compare d to 
abundances observed in diffuse clouds dLiszt & Lucasl2002l e.g. 
X{Cli) x 4.3 10 in ^ Per). On the other hand, the relatively 
constant abundance ~ 1.0 10"** at the offsets Ax < +8' is likely 
to correspond to the steady-state value in a dense core (aver- 
aged over a 4' beam). It is, however, somewhat hig her than the 
abund ances predicted by the chemistry model of iRuffle et"al] 
(Il997h where the interaction with grains is included. Then again 
co mparing our calcula ted abundances to the model presented 
bv lFlower et al.l (Il994h . we see that our abundances are slightly 
lower than the values presented in the model results where the 
depletion of oxygen 6o - I for a cloud with high visual optical 
depth. Our results would most likely match Flower's model if a 
slightly higher oxygen depletion factor is used as a model where 
(5o = 10 produces much lower abundances than our measure- 
ments. 

4.3. The structure of the filament in the light of the SCUBA 
map 

In addition to providing us with an estimate of the hydrogen 
column density in TMC-1, the SCUBA maps can be analyzed 
from the perspective of examining the physical dimensions of 



Table 2: Properties and closest correspondence with Hirahara 
cores of the four clumps found using the clumpfind algorithm 
on the 850 jum SCUBA intensity map (smoothed to 30" resolu- 
tion). 



Hirahara 


Mass [Mq] 


FWHM XxYl'] 


X:Y 


E [g/cm^] 


A+B 


3.9 


4.2 X 1.5 


2.8 


0.042 


C 


4.2 


3.6 X 14 


2.6 


0.045 


D 


3.4 


2.3 X 1.8 


1.3 


0.045 


E 


5.3 


3.9 X 1.9 


2.0 


0.041 


SEtip 


1.5 


2.9 X 14 


2.2 


0.032 



the clumpy structures seen in them. That is what we intend to do 
in this subsection. 

The clumpfind algorithm (IWilUams et al.|[T994l) was used to 
identify clumps, i.e. relatively isolated density en hancements, 
in the 850 fim SCUBA map of the filament from iNutter et aH 
(l2008h . smoothed to correspond to a 30" beam. The filament 
shows rather small intensity variations. The best agreement with 
the visual inspection is obtained by setting the intensity thresh- 
old to 0.15 Jy/beam (twice the value of the rms noise in the 
neighbourhood of TMC-1) and the stepsize to 0.075 Jy/beam 
(Icr). With this setting the algorithm finds five clumps in the 
filament; the properties of these clumps are listed in Table |2] 
Th e closest c o rrespo ndence with the "Hirahara cores" from 
Hirahara et al.l (Il992h is also indicated. The int ensity stepsize 
used here is half the value recommended by IWilliams et al.l 
(|I994|).Ifa2cr step is used, the program identifies 3 clumps cor- 
reponding to the Hirahara cores A+B, C+D, and E (extending to 
the SB tip). The FWHM sizes of the clumps in x and y directions 
and the aspect ratios are calculated in a coordinate system tilted 
by 53.6° with respect to the R.A. axis 

The total mass of the filament is 18.3 Mq, and the average 
surface density is < E >= 0.042 g cm"^. Here we have assumed 
a constant dust temperature rd=10 K. The previous determina- 
tions of the the gas kinetic temperatures towards various posi- 
tions in the cloud, includin g the CP and NH3 p eaks, lie invariably 
in the vicinity of 10 K fe.g. lPratap et alJ[l997l) . and we think this 
value represents well the dust clumps seen at a 30" resolution. 
The somewhat higher temperatures derived in Sect. 13.31 by com- 
parison between the SCUBA 450/850 maps and Ay, relate to a 
228" beam which encompasses also core envelopes. The clumps 
can also be identified as five intensity maxima in Fig. [Taj 

Even though the locations of the Hirahara cores are always 
in the vicinity of a clump on the filament, the actual peaks of the 
clumps and the CCS line intensities do not overlap; both of the 
850 micron intensity peaks corresponding closest to cores C and 
D are shifted about 2' towards each other in relation to the CCS 
intensity maxima closest to the two cores. These spatial offsets 
are larger than the FWHM of the CCS observations. 

Applying the formula quoted by iHartmanii! (l2002l) for a fila- 
ment with a critical line density, one finds that the average sur- 
face density 0.042 g cm - implies a Jeans length of ~ 0.16 pc 
or 4' at a distance of 140 pc. The corresponding mass is 2.7 
Mq. The long-axes of the clumps and the avarage separations be- 
tween their centres of mass (~ 4.8') are comparabable to the de- 
rived Jeans length. This suggests that the clumps identified here 
are formed as a result of thermo-gravitational fragmentation. 
The full resolution SCUBA maps and the previous high-angular 
resolution molecular line maps of Langer et al. (19951) suggest 
that the fragmentation has proceeded to still smaller structures, 
so that each clump is likely to contain several dense cores. 
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Finally, we note the overall structure of TMC-1 does not re- 
semble the model presented in iHanawa et al ] (11994 . where the 
fragmentation of the filament is supposed to have started from 
one end (from the northwestern tip in this case). On the SCUBA 
map we see a roughly symmetrical structure, the apparent centre 
lying between cores C and D. The structure bears resemblance 
to the fin al density profiles resulting from the simulations of 
iBastien e"t al. (199]j) for the gravitational instability in an elon- 
gated isothermal cylinder 

5. Conclusions 

We report on the detection of a CH abundance gradient in TMC- 
1 with X(CH) increasing from ~ 1.0 10"^ in the northwestern 
end of the filament to ~ 2.0 10"^ in the southeastern end. The CH 
column density maximum lies close to cyanopolyyne (CP) peak 
in the southeastern part of the filament. As it belongs to the most 
simple carbon compounds, CH is likely to reflect the distribu- 
tions of C 1 or C^. Using previously published OH observations, 
we also show that the fractional abundance of OH increases 
strongly from the northwest to the southeast. The OH/CH abun- 
dance ratio lies in the range ~ 16 - 20. 

In the northwest, the CH abundance stays remarkably con- 
stant at around 1.010 The value is clearly higher than the 
steady-state values predicted for dense cores by pure gas-phase 
chemistry models, but the agreement is better with models in- 
cluding gas-grain interaction. We performed calculations using 
the UMIST Database for Astrochemistry reaction network and a 
homogenous physical model appropriate for sources such as the 
TMC-1 ridge. The best agreement with observations for the pure 
gas-phase model was found at cloud ages of several 10^ years, 
but the predicted CH and OH fractional abundances are a fac- 
tor of 3-5 lower than the observed values. Agreement is better 
with a gas-grain model at this age. The latter model reproduces 
the observed abundances at late stages where molecular deple- 
tion is significant. However, in view of the large beam used in 
the observations, and the likely encompassment of low-density 
gas, it is doubtful that the derived abundances reflect conditions 
in depleted cores. 

The dust continuum map from SCUBA dNutter et al.ll2008l 
smoothed to the 4' resolution of the CH observations) suggests 
that in the vicinity of the CP peak the filament has an exten- 
sive envelope, in contrast to the neighbourhood of the NH3 peak 
which appears to be more compressed. As the CH and OH abun- 
dances are expected to be anticorrelated with density, their larger 
abundances in the southeast are likely to be due to this low- 
density envelope. The largest fractional CH and OH abundances 
are found at the southeastern end of T MC- 1 where we are a p- 
proaching the "C 1 peak" discovered bv iMaezawa et al l (ll999h . 

The total H2 column density distribution along the TMC- 
1 filament was calculated by combining the SCUBA 850 and 
450 fim maps with an Ay map from 2MASS data. The distribu- 
tion peaks between the CP and NH3 maxima. In fact, the high- 
resolution SCUBA map at 850 yum shows a rather symmetric 
structure with the apparent centre between the Hirahara cores C 
and D. The sizes and typical separations of clumps correspond 
to the Jeans' length in an isothermal filament with the observed 
surface density. The main difference between the two ends of the 
filament seems to be the presence of the low-density envelope in 
the southeast, which probably is an indication of an early stage 
of dynamical evolution. As lower densities result in slower reac- 
tion rates, this finding conforms with the suggestion of several 
previous studies that the southeastern end of the cloud is also 
chemically less evolved than the northeastern end. 
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